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We present a theoretical analysis of the properties of low-dimensional quantum antiferromagnets 
in applied magnetic fields. In a nonlinear a model description, we use a spin stiffness analysis, a 
1/N expansion, and a renormalization group approach to describe the broken-symmetry regimes of 
finite magnetization, and, in cases of most interest, a low-field regime where symmetry is restored by 
quantum fluctuations. We compute the magnetization, critical fields, spin correlation functions, and 
decay exponents accessible by nuclear magnetic resonance experiments. The model is relevant to 
many systems exhibiting Haldane physics, and provides good agreement with data for the two-chain 
spin ladder compound CuHpCl. 



PACS numbers: 75.10.Jm, 75.30. Cr, 75.40. Cx 
I. INTRODUCTION 

The importance of low-dimensional spin systems in re- 
vealing fundamental quantum mechanical properties has 
been recognized since Haldane's conjectured concerning 
the effects of quantum fluctuations in integral- and half- 
integral-spin antiferromagnetic (AF) chains. The key 
quantity is the topological term arising from the quan- 
tum spin phase, and similar considerations have since 
been extended to anisotropic and higher-order interac- 
tions, higher dimensions (planes), and coupled chains 
(spin ladders), largely driven by materials and experi- 
ment. Magnetic flux within a sample can alter the effec- 
tive spin magnitude, and features visible in an applied 
field include spin gaps and magnetization plateaus.tl 

Recent progress in metalloorganic synthesis has fur- 
nished new classes of low-dimensional antiferromagnets 
(AF) whose small exchange constants make their full 
magnetization response, including saturation, accessible 
to laboratory fields. The measured magnetization curves 
of some such materials, thought to be prototypical quan- 
tum magnets, in fact exhibit certain features which are 
rather classical. Examples on which experimental intefc 
est has focused are the Haldane (5 = 1) chain NENP,lI 
the planar system CFTDp and the two-chain spin lad- 
der CuHpCl.tffl We will focus on the ladder geometry of 
CuHpCl, the best-characterized sample in recent litera- 
ture. 

We consider the quantum AF system in an exter- 
nal magnetic field using the nonlinear a model (NLsM). 
While this widely-applied treatment is in fact semiclassi- 
cal, being truly valid only in the limit of large on-site spin 
S, it has in the past formed the basis for many funda- 
mental deductions concerning the quantum limit of an- 
tiferromagnetic spin systemsoij. We will demonstrate 
here its validity in the case of all effectively integral-spin 
quantum systems in appreciable magnetic fields, and pro- 
vide justification for this result in terms of suppression of 
quantum fluctuations by the field. Such systems display 
an ideal quantum phase transition, driven by the applied 
field, between a disordered regime with gapped spin exci- 



tions at low field and a quasi-long-range ordered regime 
with gapless excitations at higher field. 

This model presents an ideal example of symmetry 
breaking in condensed matter systems. The 0(3) symme- 
try of the spin Hamiltonian, and of its low-energy descrip- 
tion in terms of the NLsM, is broken on a purely classical 
level. However, at low fields this symmetry is completely 
restored by strong quantum fluctuations, a well-known 
property of the NLsM, which has been used to illus- 
trate the effect pf asymptotic freedom (or "confinement 
of excitations" ). 113 In contrast, quantum fluctuations may 
be suppressed by fields which are sufficiently strong but 
remain experimentally accessible, with the consequence 
that the full 0(3) symmetry cannot be restored, and the 
spin system is reduced to the lower XY or 0(2) symme- 
try. This XY symmetry, in the plane perpendicular to 
the applied field direction, remains sufficient that there 
can be no true long range order, but instead a quasi- 
long-range order characterized by spin-spin correlation 
functions which decay in space-time with a power law 
form. 

The outline of this paper is as follows. In Sec. II we de- 
rive the form of the NLsM in a magnetic field, using the 
geometry of the two-chain ladder system. With a view to 
experimental comparison we include magnetic saturation 
by applying a total-spin constraint. In Sec. Ill we analyze 
the spin stiffness of the model, which allows us to deduce 
the general behavior, correlation length, spin gap and 
critical field at zero temperature, and also the effects of fi- 
nite temperature and system size. In Sec. IV we consider 
the 1/N expansion, which is well suited for describing the 
low-field, disordered regime. Sec. V contains an exten- 
sive renormalization group (RG) study of the model over 
the entire field range. We derive the coupled RG equa- 
tions for the 1+1-dimensional (1+ld) system and present 
their general solution. This allows us to show how the 
equations return the physics of symmetry-breaking and 
-restoration in the previous paragraph, and to discuss 
the issue of renormalization and gauge-invariance. We 
compute in Sec. VI the magnetization of the model us- 
ing a high-field expansion, and in Sec. VII calculate spin 
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FIG. 1. Representation of two-chain ladder system and co- 
rdinate axes. J = J^ is th< 
interchain exchange coupling. 



ordinate axes. J — J x is the intrachain and J' = 2Jj, the 



correlation functions for the low- and high- field regimes. 
In Sec. VIII we compare our results with experimental 
data from magnetization and nuclear magnetic resonance 
(NMR) spin-relaxation measurements. Sec. IX summa- 
rizes our conclusions, and discusses the variety of systems 
and physical problems to which the formalism is applica- 
ble. 



II. NONLINEAR a MODEL 

We consider the ladder geometry of the compound 
CuHpClJj The Hamiltonian for the system in a magnetic 
field b = 5/isB may be written as 

H = 2_^ [J^rn,i ' S„j : i+i + J Sl,i • S2,i + b • S m ,j) , 
i;m— 1,2 

(2.1) 

where J is the superexchange interaction between spins 
in each chain and J' the interchain, or ladder "rung" , 
interaction. We choose iV^-site chains with periodic 
boundary conditions in both directions, and the geom- 
etry shown in Fig. 1, such that J x = J, J y = J'/2, 

B = (0,0,5), and S miWx +i = S m ,*. Followina.ti^rSjtan- 
dard procedure in the derivation of the NLsM,EiEirL3 the 
coherent-state representation of the spin S„ M ; is written 



as o & t r 



s[(-i) 



1+771 



n„ 



}lm,i], where n m ^ is a stag- 



gered spin (the Neel vector) with unit magnitude, and 
lm,i describes the spin fluctuations perpendicular to n. m ,i- 
On proceeding to the continuum limit, the Hamiltonian 
is expressed in coherent-state representation as 



1 



H = dr \o s2 E J ^ ^ n ) 
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(2.2) 



The full Euclidean action Se = Swz+f drH in space 
r = (x,y) and inverse temperature t (/3 — l/ksT) con- 
tains a Berry-phase, or Wess-Zumino, term whose origin 
lies in the solid angle subtended by the fluctuating quan- 
tum spin phases: 



S wz = i— drdr 4> - a\ ■ (n A ri) + AiriS (Pi + P 2 ) . 

(2.3) 

Here, <fi is the azimuthal angle, and P\ = P2 — 
(1/4-7t) J drdx (n A h) • d x n are topological terms which 
give the Pontryagin index on each chain whcp-^n = 0, 
as is the cas e in a ladder of only two chains. IB The last 
term in Eq. J2.3| ) is therefore z(47rPi)25, demonstrating 
that the system will have integral-spin characteristics for 
any value-pcf S, and the topological term may thus be 
ignoredBO 

Integrating out the fluctuations 1 about the staggered 
spin configuration of n, subject to the orthogonality con- 
straint n • 1 = 0, yields the classical solutionolifl 
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{i(nAn) + [b-n(b-n)]}, 



(2.4) 



where J — J x + J y . Finally, the corresponding action for 
the quasi-one-dimensional (ID) ladder system, in 1+1 
Euclidean dimensions denoted by /1, is 



Se = Tr- I dr dx 



{(cVn) 2 -[b 2 -(n.br 



2ib-(nAh)|, (2.5) 



where g = (2/N y S)\/j/J x is the bare coupling constant, 
and the integral over r is from zero to Lt — C/3, with 
c = 2Sa\f J X J the effective spin-wave velocity. (We use 
h = 1 throughout.) Summation over fi is implied in the 
square. We have left explicit the number of chains N y 
in the ladder. For CuHpCl we have that J y > J x and 
N y = 2, whence the rung coupling may be treated as 
effectively rigid. T he fo rm of the NLsM in an exten 



field given by Eq. ( |2.5|) has been derived previously, LLI 
and its implications considered in the low-field limit x3 
We will examine in this work its properties at all field 
strengths. 



The action (2.5) may be cast in a physically more 
transparent form by expressing the staggered spin in 
spherical coordinates, n = (sin 8 cos (f>, sin 9 sin (j>, cos 8). 
From the resulting Lagrangean density 
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+ {d,8f\ (2.6) 



one observes that the field induces a hard-axis anisotropy 
through the term (n • b) 2 , which makes the spins tend 
to align with the Neel vector n in the plane normal to 
b. For strong fields, deviation of n from this plane is 
small, and expansion to quadratic order in •d = 8 — n/2 
is valid. Taking the field along the z-axis, we find that 
the Lagrangean 
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di + (j> + ib)' 



(2.7) 



separates into contributions from in-plane and out-of- 
plane fluctuations in the high-field regime, \b\ > p = u. 
We will use Eq. (|j), and return to the condition for its 
validity, below. 

Finally, for a real magnetic system with small spin 
quantum number S, application of a sufficiently large 
magnetic field will cause total spin alignment, or satu- 
ration, with a maximum magnetization given by M s = 
g/j,BSN s , where N s (= N X N V here) is the total number 
of sites in th e sy stem. This effect is not contained in the 
NLsM (Eq. (|2.5|)) discussed above, where S is required 
to be large (semi-classical). Furthermore a saturation of 
the spin would correspond to large |1| and sma ll |n |, ren- 
dering invalid the expansion of S (below Eq. (2.1)). We 
will include saturation by placing a limit, corresponding 
to the desired experimental spin magnitude, on the spin 
component (S z ) projected along the field. By this tech- 
nique, the full spin vector remains primarily in the plane, 
justifying the treatment of the uniform spin component 
|1| as a small variable. 

The succeeding experimental comparison (Sec. VIII) 
will be simplified by treating M a as a simple cutoff, 
but we describe here briefly this systematic inclusion 
of saturation in the model by the application of a con- 
straint on the total spin. In addition to the constraints 
n • 1 = on the uniform and staggered spin compo- 
nents, and n 2 = 1 on the magnitude of the staggered 
spin, we introduce a further constraint enforcing conser- 
vation of spin. For a field applied along the z axis, the 
conserved quantity is the z-component of the total spin 
S z = {{fl}\S z \{fl}) -* Sfdxl z , where the continuum 
limit is taken in the final step. While the previous con- 
straints are purely local, in that they must be enforced at 
each spacetime point, the spin constraint is local in time 
but global in space. All of the constraints are treated by 
inserting (functional) delta functions into the measure of 
the propagator, in the form 

<{0} 6 |e-^|{0} Q ) 

= fvn(x,T) fvi(x,T)e- fc[n - l] e- SE[n - l] , (2.8) 

where the function /c [n, 1] specifies the constraints using 
Lagrange multipliers. Integrating over the fluctuating 
spin field 1, which is now constrained, and over all the 
Lagrange multipliers, one obtains the Euclidean action 

<S[n] = S W z + — drdx (d^n) 2 

_ylmw£_ 81l ^ (2 „ 

where D[n] = Jdx (n A n) z and A[n] = Jdx (1 — n 2 ), and 
we have introduced the scaled spin o = S Z /(SN S ), —1 < 
a < 1. In this action the field b is coupled only linearly to 
the spin a, as in ferromagnetic systems, and in contrast 
to Eq. (|]^), where it appears quadratically. This linear 
coupling is a direct consequence of the spin constraint, 
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FIG. 2. Ladder system with twisted staggered moment. 



and one may recover the previous action (Eq. (2.5)) by 
allowing a to fluctuate in time, relaxing the condition 
\a\ < 1, and integrating over ct(t). Evaluation of the 
partition function with s atu ration, requires summing the 
constrained action (Eq. J2.9| )) over all degrees of freedom, 
which include the allowed values of \S Z \ < N S S. This 
procedure limits M to M s at large fields, and provides 
a consistent description of a second phase transition, to 
saturation, which occurs at high fields. In Sees. III-V 
we will concentrate on the symmetry-breaking quantum 
phase transition (Sec. I), which occurs at a lower field 
on the order of the zero-field gap (below), and so the 
discussion there will not involve further consideration of 
the spin constraint. 



III. SPIN STIFFNESS 

To gain initial insight into the effect of the magnetic 
field, we consider the spin stiffness of the ladder system 
using the method of Ref. [Eg|£3 We take (Fig. 2) the stag- 
gered spin configuration to be subject to a twist tp in the 
plane, n(r, 0) = (1,0,0) and n(r, L) — {cos ip, simp, 0). 
This is conveniently represented by n = "R,{x)(T, where 
1Z is a matrix for rotation by ipx/L about the z axis, and 
er = ct(t,x) a 3-component vector with value (1,0,0) at 
the spatial endpoints. The cross-term in (<9 M n) 2 and the 
third term in Eq. (2.5) are total derivatives in x and r 
respectively, and will be discarded for the present pur- 
poses. This is equivalent to considering only paths with 
no winding in Euclidean space, and gives the action the 
form 



Sl=~ l<lr<l-r^,) /l( Tr- + [±2-^)0-ni) 
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(3.1) 



To 1-loop order in g, the spin stiffness is given by 

As 
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(3.2) 



where p° = c/2g is the classical (bare) value, and the sum 
includes both quantum and thermal (through the finite 
"length" Lt) corrections to first order in g. 

We consider first the low-temperature, or "quantum" 
regime Lt S> L. We evaluate the summation between the 
spatial limits tt/L and ir/a (this change to open boundary 
conditions does not affect the result), and introduce the 
"magnetic length" 



L m = nc/b 
to obtain the expression 



(3.3) 
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In 



a/L + y/(a/£) 2 + (a/L m y 
4 (l + y/l + (a/L m ) 2 ) 



(3.4) 



The system length L in Eq. (3.4) may be substituted by 
a correlation length £, beyond which segments of the lad- 
der behave independently, and £ then computed from the 
condition p s = 0. We emphasize that p s is to be consid- 
ered as a local stiffness, meaningful only on length scales 
L < £. In the zero-field limit (L m — ► oo) one obtains the 
result £o = Aae 27T / g = Aae a7rS , where a — \/,J x /J intrp-. 
duces a dependence on the ladder coupling constants, tZl 
and A is a nonuniversal constant of order unity which 
depends on the integration cutoff. The general solution 
can be expressed as 



£(£) = 2L* m /[l -(L* m /L m ) 2 ] 
= £ /[l-(B/B*)*], 



(3.5) 



from which we see that the correlation length £,(B) di- 
verges at a critical field B* , where L* m = asinh(27r/g). 
For fields B < B* , the system has only short- 
range correlations, and the finite correlation length 
may be written as £(B) — ae Q7r5 , where S — 
S [l - (g/2n) In (l - (L* n /L m ) 2 )] is a growing value of 
the effective spin. There is no spontaneous breaking of 
the 0(3) spin symmetry. For B > B* , the field etu 
forces a long-ranged correlation throughout the system,Eil 
and it is most convenient to write the spin stiffness as 
Ps = p° [1 — (<?/27r) ln(l + L m /a)], a quantity which re- 
covers the bare value p° s as B — ► oo. Finally, at the 
transition B* , it is clear from Eq. ( p^g ) that the diverg- 
ing correlation length corresponds to the closing of a gap 
A ~ £ _1 to spin excitations with the form 



Aoc 1- (B/B*Y 



(3.6) 



This situation is summarized in Fig. 3. While A is the 
gap to all fluctuations in the low-field regime where 0(3) 
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FIG. 3. Schematic behavior of spin stiffness and correlation 
length, or spin gap, with applied field. 

spin symmetry is maintained, we note that at high fields 
not all the modes are gapless. In this regime (strictly 
B ^> B*) the symmetry is l owe red to 0(2), and while in- 
plane ((/>) fluctuations (Eq. (2/7)) are massless, excitation 
modes (§) in the direction of the field, or out of the plane 
it enforces, are gapped with "mass" b. 

We may further employ the spin stiffness as a means 
of characterizing the behavior of a system of finite size 
L, and at finite temperatures, for different values of the 
magnetic field. This is essentially a matter of comparing 
the corresponding length scales L, Lt and L m E£l. We 
present results for the cases of i) weak and ii) strong fields 
in the following three regimes: a) quantum, i) L m , Lt 3> 
L and ii) Lt » L» L m , b) classical renormalized, i) 
L mi L ~^> Lt and ii) L ^> Lt 3> L m , and c) classical, i) 
L m , L ^> a ^> Lt and ii) L ^> L m ,a ^> Lt- Following 
Rcf. [19] we quote the stiffness in dimensionless units of 
p s jc — (\/2g){p s / p° s ). For weak fields we find 



c 



(3.7) 

£0 denotes the quantum correlation length defined above, 
£ c l = 6Lt/ ' g is the classical correlation length, and 
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(c^i + ^T+(L T /L m y 



(3.8) 



with C ~ O(l) another nonuniversal constant, contains 
the renormalizing effect of quantum fluctuations on £ c i. 
Thesfi-jesults are all directly analogous to the zero-field 
case,c2l with the field entering as a quadratic correction 
where relevant. In a strong applied field (ii), 
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(£t(B) - L m )/12LT, 

(Cci - L m )/12L T , 



L T > L > L m 



(3.9) 



where now the renormalized classical correlation length 
has the field-dominated form 
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Id ( Cf- V Vl + (L m /L T ) 2 + L m /L T 



(3.10) 



The general characteristic of all the above results is a log- 
arithmic dependence of the spin stiffness on the length 
scales of system size, temperature, and field where quan- 
tum fluctuations are important, turning to a power-law 
dependence in the classical regime. In the high-field sit- 
uation the scale L m of the field repla ces Lt as the de- 
termining parameter (cf. Eqs. (3.7, 3.9)). Equivalent 
results for finite systems may also be deduced from the 
1/N expansion, which we consider in the following sec- 
tion. 

We will later compute the magnetization M = 
dF/dB\^, = o, both in the thermodynamic limit and for 
systems o f fin ite size, from the same type of free energy 
as in Eq. (3.2). This undertaking is simplified in the light 
of insight gained from 1/N and RG analyses, and is thus 
deferred to Sec. VI. At this point we note only that the 
magnetization calculated in this framework is not zero 
in the low-field regime, as would be required of a sys- 
tem with a spin gap and no broken symmetry. This is a 
consequence of the approximation, implic it in Eq. (2.7) 
and in deriving the free energy F(b) (Eq. ( ft.2| )), that the 
field b be large on the scale of \<j>\. We may neverthe- 
less conclude that analysis of the spin stiffness itself is 
qualitatively revealing of many properties of the system. 



IV. 1/N EXPANSION 

To address the situation in tlus weak-field regime, we 
adopt instead a 1/N expansion)!!] which is expected to be 
appropriate in describing spin-gap phases. Here the stag- 
gered spin n is assumed to exist in an TV-dimensional spin 
space, in which only the direction n z is selected by the 
magnetic field, and a controlled expansion may in princi- 
ple be performed in the small quantity 1/N to compute 
fluctuation corrections to the saddle-point solution. The 
relevant parts of the action are (cf. Eq. (3.1)) 



S t 



2.9 



dr dx 



(d,n) 2 -b 2 (l 



-iAn'-l 



(4.1) 



where b denotes b/c, and the constraint that the spin n 
have unit magnitude is made explicit with the Lagrange 



multiplier iX. Integration over n and taking the func- 
tional derivative of the resultant effective action with re- 
spect to A yield an equation for the saddle-point value of 
the Lagrange multiplier, 



3 L- 1, 



1 



k 2 + iX + b 2 '' 



(4.2) 



in which the J3-field term is found to appear only at 
0(1/N). iX appears as a mass, or cutoff term in mo- 
mentum integrations, and is thus an upper length scale 
for cooperative processes in the system, or simply a cor- 
relation length (inverse excitation gap) . Thus by writing 
the saddle-point solution as iX — c 2 n 2 /£ i (B) 2 , and car- 
rying out the summation at low T from to ir/a, this 
correlation length is determined from 
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y/(a/Q 2 + (a^L-y 
1 + v/1 + (a/0 2 + (a/L m y 



(4.3) 



whic h sh ould be compared with p s — emerging from 
Eq. Q . 

In the weak-field regime, £ <C L m , one finds 



£ ~ (l/2)ae 27T/Ng [l + (a 2 /8NL 2 Je 4n/g ] 



(4.4) 



This is an 0(1 /N) correction to the previous result, 
but there is also a difference of a power of 1/N in 
the exponent. Taking the field derivative, the result 
<9b(c7r/£) 2 = —2b/N ensures that the magnetization con- 
tribution from the k summation terms in F are identi- 
cally zero to 0(1/N). At the same order, the field deriva- 
tive of the k-independent term b 2 +c 2 n 2 /^ 2 is 2b(l—l/N). 
The behavior required of a gapped system is obtained in 
the weak-field regime on making the well-known identi- 
fication, deduced by comparison with exact RG results 
(see Ref. |l(] and below) N — > N — 2. Returning to the 
physical situation of the 0(3) magnet, the magnetization 
is zer o at all fields B < B* . The saddle-point solutio n fo r 
£(.B) ( |4.3D becomes precisely that deduced from Eq. (3.4) 
(up to a negligible term retained from a lower cutoff) , 
with the same field B* for divergence of £. For the 0(3) 
case we thus deduce the expected gapped state, with no 
breaking of spin rotational symmetry. 

For completeness, we note the strong-field results 
of the 1/N treatment. When the magnetic length 
is small (L m <C £); the correlation length £ ~ 
(l/2)ae 27T / 9(N ~ 1 '>(a/L m ) 1 /( N -V has a direct field depen- 
dence. The magnetization is given to leading order by 
M ~ -C + 2bL 2 m /[(N - 1)£ 2 ]. The quantity appearing 
in these expressions is not N but N — 1, indicating again 
the lowering of symmetry in the field direction above a 
finite critical field B*. However, these high-field results 
are not meaningful in the 0(3) case. 



In summary, the spin stiffness and 1/N techniques 
taken together yield a consistent picture at both weak 
and strong magnetic fields. However, while the Ansatz 
N — > N — 2 has been applied previously]!!] and leads to 
consistency between the results of this and the previous 
section, it cannot be taken to be entirely satisfactory. For 
this reason we pro ceed by applying the RG technique to 
the model of Eq. (§j|. 



V. RENORMALIZATION GROUP 

We consider next a renormalization-group (RG) study 
of the NLsM in an applied field. This approach yields 
meaningful results over the full parameter range. We 
adopt a standard Wilson momentum-shell treatment, 
further details of which are presented in App. 1. In brief, 
the action in the form 



2g 



<h 
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dx [{dfi\ 



i) 2 - 6 2 (1 - n\) + Aibn x hy] 



(5.1) 



may be reexpressed in terms of variables <f>, representing 
in-plane fluctuations of n, and ^fgo z — n Zl representing 
out-of-plane fluctuations. The latter is chosen to facili- 
tate a perturbative expansion in g, which yields 



C E = 4-(A>-b 2 )-±a z (-dl + b 2 -A>)a z 



0(g). 

(5.2) 



A' = (d^cj)) 2 + 2ib(j) denotes in-plane fluctuation terms, 
which, because <f>(r, x) is assumed to vary slowly, can be 
taken to be a small constant (no fast Fourier modes) in 
the momentum shell 7A < |k| < A (7 — > 1) around the 
finite upper cutoff A = n/a set by the l attic e length scale. 
By expanding in A' , the form of Eq. ( |5.2j ) is recovered, 
but with new coefficients g{a') and b(a') z given at 1-loop 
order by 



-L-tx—L 

9o -dl H 
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9(a') 9o " -d 2 +b 2 ' 
b{a') 2 = bl- ff0 Tr' '\n(-dl + bl) 



The integrals represented by the partial traces (Tr') are 
performed over an isotropic momentum shell in 1+ld, 
with At = A. Alterations arising from considering in- 
stead a rectangular momentum shell, with separate inte- 
gration over space and inverse temperature, would affect 
only non-universal prefactors. 



Inspection of the equation (5.3) for rcnormalization of 
the coupling constant g indicates already the effect of the 
field. Performing the momentum-shell integral yields 
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(5.5) 



for a < a' < £ = min[X,Lx], assuming (see below) 
that the flow of L m is at most weak. For weak fields 
L m ^ £ we see that go/g(C) = 1 — (go/2n) ln(£/a), 
i.e. g is renormalized by the full effects of quantum 
fluctuations. For strong fields C » L m ^S> a, we have 
go/g(C) = 1 — (go/2ir) ln(L m /a), and the coupling con- 
stant g is independent of the system length scale C In 
this case renormalization is weak, and for very strong 
fields L m <C a it vanishes (go/g(C) = 1), demonstrating 
the suppression of quantum fluctuation effects by high 
fields. 

Following the standard RG procedure, we define the 
flow parameter I = ln(a'/a), whence dl = da/ a. Tak- 
ing the derivative with respect to I of Eqs. (5^,5^4), we 
obtain the differential form of the coupled RG equations 
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These represent an extension of the usual RG equations 
to include a mag netic field B, contained in (3 = a'b(a') (= 
a'b(a')/c). Eq. ( |sli| ) is the conventional "/3-equation" , 
the terminology given to the renormalization of the cou- 
pling constant, but with an additional field term in the 
denominator. Clearly, a strong field restricts flow to the 
strong-coupling (disordered) limit, and, as observed in 
the previous paragraph, acts to suppress quantum fluc- 
tuation effects. This points to a "deconfinement of ex- 
citations" at suitably high field. Eq. ([5J) expresses the 
renormalization of the magnetic field term b — gfJ-sB 
with dynamical exponent z — 1 (from the first term, 
which arises purely from the presence of a' in the quan- 
tity /3), but with additional, logarithmic suppression of 
this flow at strong field and coupling (second term). 

At this point some commentary is in order regarding 
the above RG equations. The NLsM in the presence of 
a magnetic field may be expressed in the gauge theoretic 
forn£3t 
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(5.8) 

with the field term cjijlbB/c in the tale of a gauge poten- 
tial. From this it has been argueda that the field term 
should scale only with t emp erature (the z = 1 behavior 
of the first term in Eq. (5/7)), and cannot be renormal- 
ized separately, for example by interactions, as this would 
violate the gauge invariance. This sca ling hypothesis 
stands in contrast to the result (Eq. (5.7)) of the Wilson 
momentum-shell RG approach taken above, where the 
second term indeed gives a field renormalization, albeit 
one which is logarithmically weak. However, the present 
situation is not unique: similar breakdown of this scal- 
ing occurs in other cases, particular examples being the 
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FIG. 4. RG flow diagram for g and b. Strong- and 
weak-coupling regimes are separated by separatrix s. 



Fermi liquid and the Luttinger liquid, also considered in 
Ref. B3, where there is a temperature-independent renor- 
malization of the field termO'El We note also that on 
taking into account magnetic saturation (Sec. II) in the 
present model, additional field terms would arise which 
preclude the use of a gauge-invariance argument. Finally, 
the qualitative features of the discussion are not affected 
by the logarithmic term in Eq. (5.7), because its field 
renormalization is weak. 



In analyzing the content of the RG equations (5.6,5.7), 
we concentrate first on the fixed points, in order to obtain 
a qualitative picture of the RG flow diagram. Seeking a 
fixed point by a weak-field expansion around 0* — 0, we 
find 
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which may be solved to yield 
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The fixed point (<?* ,/?*) = (oo, 0) is evidently stable if the 
flow is stopped at I* — 2*k j go- The system will flow to 
this strong-coupling regime if the starting value (3~o is suf- 



ficiently small. The length scale £* 



J» 



at which the 



flow stops may be compared with the spatial and ther- 
mal dimensions (L, Lt) of the system to calculate directly 
the effects of finite size and temperature (cf. Sec. III). 
At strong fields (/3 — ► oo), dg/dl = 0, or g = <?o,_and sim- 
ilarly d\n(3 2 — 2dl, from which it follows that b(l) = bo, 
i.e. in this regime neither the coupling nor the field is 
renormalized. 



The RG equations (5.6 5.7) may be solved numerically 
in the form 
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where x = 2tt / g is the inverse coupling constant and 
y = ab corresponds to the field variable b/c only. Solu- 
tion leads to the flow diagram in Fig. 4, which has the 
following interpretation. The regime (i) of weak initial 
B-field is the strong-coupling phase, with confinement of 
(gapped) excitations. Here, the assumption (underlying 
the perturbative RG treatment) of small g becomes in- 
consistent, but one may still deduce the critical length 
scale £* , and that the magnetization M = 0. In this re- 
gion 0(3)-symmetry is restored by quantum fluctuations, 
which may thereby be considered as "screening" the mag- 
netic field. In contrast, the regime (ii) of strong initial 
B-field corresponds to weak coupling, where g and b are 
only weakly renormalized. Here, the excitations are de- 
confined on a length scale £(B) whose flow is governed by 
B. In this region, quantum fluctuations are suppressed 
by the magnetic field, and the broken 0(3) symmetry 
cannot be restored. 

The properties of the broken- and unbroken-symmetry 
phases may be further contrasted by considering the cor- 
relation length £ in each regime. This is a physical quan- 
tity and so does not change under the RG flow, whence 
d^/da' = and 
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d(3 2 d£ 
dl dp 2 
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At low fields, we find that — <9£/£ ~ 2irdg/g, which has 
the solution £o = ^e 2 ^ 1 / 90-1 / 9 '. If one assumes that 
£ — ► a (the bare lattice constant), as the system flows to 
the strong-coupling limit (I /go — * 0), then £o = ae 2lx l 90 
returns the same, finite physical correlation length as in 
the previous analyses. For high values of the initial field 
we have — 2<9£/£ ~ d/3 2 //3 2 , which leads to the relation 
£o = ££/a (bo is invariant) under the RG flow. Thus, for 
any finite £ during the flow, the bare correlation length 
£o is the system size L or Lt, corresponding to the quasi- 
long-range order expected of an XY system in one spatial 
dimension. 

Before concluding our RG analysis, we also consider 
its qualitative predictions for the system magnetization, 
a measurable quantity. This may be calculated as 
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Once again, 
and dfl/dfto 



for weak fields we find that dg/d(3o 



o 



& \J \. — gol/2n, where the latter expres- 
sion will introduce the characteristic length scale C of the 
short-range correlated phase. Now for any F(b) analytic 
in 6, 
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from which it is clear that the magnetization M(Bo) van- 
ishes in the scaling limit C — ► C* corresponding to strong 
coupling (i.e. for all bo sufficiently close to the fixed-point 



value bo = 0. Within linear response, this vanishing of 
the magnetization M = x±B corresponds to the vanish- 
ing of the susceptibility x± ~ m (£/-k) as thaiength of the 
system increases into the disordered stateila Hence the 
RG framework returns a zero magnetization in the regime 
where quantum fluctuations ensure no spontaneous sym- 
metry breaking. At high fields the magnetization is fi- 
nite, but its saturation in the physical system limits the 
useful information which can be deduced from the above 
approach without inclusion of the saturation constraint 
(Sec. II). 



VI. MAGNETIZATION 

We comp ute the magnetization from the Lagrangean 
£e of Eq. (2.7), which has separate contributions from 
in- and out-of-plane spin fluctuations. The calculation 
is valid in the regime |6| > </> = u> quoted in Sec. II, 
and is thus effectively a high-field expansion, us denotes 
the energy of excitations of the system, which may be 
taken to be spin waves (w = ck) in the gapless regime. 
The magnetization is given by M — g/indF/db, with 
F = —j3~ x In Z obtained from the partition function Z = 
J D-dDif) ex p(— i Se). Within the approximations specified 
below Eq. (2/7), the separable nature of £e gives the 
result M — M +Mi. The out-of-plane contribution M is 



purely dynamical, in the sense that ■& fluctuations appear 
only in the form (l/2g)d G~ lr d. We may write G _1 = 
Gq 1 — X , with the bare propagator given by G^ 1 — b 2 — 
d 2 — d 2 and the corrections by X — cb 2 + 2ib<j). In the high- 
field regime \4>/b\ is small, so that contributions from X 
may be neglected. The appropriate piece of the partition 
function is readily evaluated as 
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The assumption of an excitation spectrum linear up to 
the zone boundary is easily altered, for example to a sinu- 
soidal form u>k = | sin ka\ within the Brillouin zone, which 
gives rise only to small numerical corrections. Working 
consistently to lowest order in the small parameter c/b, 
and in limit of low T, we find that this contribution is 
simply a constant, M = hN x g[iB, which corresponds to 
a spatially uniform state. Small corrections of the form 
M' Q ~ B In B arise for either type of dispersion. 

Treatment of the in-plane variables <f> at finite temper- 
atures is more involved. We first represent the relevant 
part of the Lagrangean as 
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to separate the term quadratic in b which gives the classi- 
cal, linear magnetization. A full treatment of <fi £ [0, 2ir] 



involves decompactification to the interval — oo < (f> < oo 
and restoration of periodicity by summation of the par- 
tition function over all winding number sectors m G Z. 
The phase is represented by <p(x, r) = <fi(x, r) + 27rmr//3, 
with cj> a periodic function of the inverse temperature 
((/)(x,t + f3) = <p(x,r)). This form takes into account 
the possibility of windings in the space of r, whereas 
windings in x may be safely neglected for a large system 
(L 3> Lt). On expanding the action, noting that all 

terms L drip vanish (integral of total derivative), and 
discarding all terms independent of b for the purposes of 
computing the magnetization, we obtain the correspond- 
ing partition function 
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Here 



is 



Co 
2 )/2e, 



a constant and a = N s b/AJ — N x ln(e + 
\/l + e 2 )/2e, wit h e = ire/ab, arises from the 0-linear 
terms in Eq. (|6J) . Because e is a small parameter in the 
approximation we employ, the second term in a simplifies 
to N x /2; this is also the origi n of the results for M and 
M' a . The sum over m in Eq. ( |6.3| ) is simplified on recog- 
nizing that Z may be represented by the 83 function,Lj 
in terms of which the (logarithmic) 6-derivative yields 
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The latter term is a sawtooth form, shown as the quan- 
tity Mf in Fig. 5(a), which is periodically extended from 
the interval —1/2 < a < 1/2 to all of a. When super- 
imposed on the linear part contained in the first term 



of Eq. (3.4), Mf gives the step-like form shown in Fig. 



5(b) . The width of the steps scales as 1/N S , so we see that 
Mf is the finite-size correction to the linear magnetiza- 
tion. Such effects are of considerable interest in molecular 
magnets, where th^ small number of atoms gives rise to 
well-defined steps.EJ However, Mf becomes indiscernible 
as the system size is taken to the thermodynamic limit, 
which is the condition we investigate below. In closing 
this section, we note that evaluation of the magnetization 
in the NLsM for weak fields b, and particularly around 
B* , remains an open problem. Experimental comparison 
for this regime (Sec. VIII) is facilitated by the knowledge 
(Sees. III-V) that M = below B*. 



VII. SPIN CORRELATION FUNCTIONS 

Correlation functions describing the nature of the spin 
order may be calculated directly in the NLsM framework 
from 




similar to those of a Luttinger-liauid description of the 
gapless, broke n-sy mmetry phase£3 The correlation ex- 
ponent in Eq. (7J5) has a direct dependence on the bare 
coupling constant, and may be identified as an effective 
Luttingcr-liquid exponent 



Mi 
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FIG. 5. Schematic representation of finite-size magnetiza- 
tion contribution Mi (a) to linear magnetization term (b) at 
high magnetic fields. 

(f r ^(x,r)4-(0,0)> = (-l) x S 2 (n l (x,T)n J (0,Q)) 

+ (aS) 2 (h(x,T)l,(0,0)) : (7.1) 

and used to compare the properties of the system at low 
and high fields. T T is the ordering operator in imagi- 
nary time. In the low-field regime where full symmetry 
is retained, the gap A to all spin excitations results in 
exponentially decaying spin correlations, reflecting the 
short-range nature of the quantum disordered phase. 

Above the critical field, as discussed previously, the 
gap closes and the low-energy features are described by 
an XY model for the plane normal to the field. The cor- 
relation functions in this regime are discussed in detail 
in App. 2. Because the quantities (n z (x, T)n z (0, 0)) and 
{l+(x, t)1-(0, 0)) (l± denotes l x ± il y ) depend directly 
on the gapped, out-of-plane fluctuations in $, these re- 
main short-ranged. However, (l z (x, t)I z (0, 0)) is propor- 
tional to d 2 ((j)(x,T)<j)(0, 0)), which leads to a power-law 
(or "quasi-long-ranged") decay of the correlation func- 
tion parallel to the field direction, 
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We have omitted here a constant term proportional 
to the square of the local magnetization. Similarly, 
(n + (a;,r)n_(0, 0)) = (exp(i</>(x, t) — i</>(0, 0))), which 
gives 
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a separate power-law form for the correlation function 
perpendicular to the field direction. B and C are con- 
stants. The results for the long-ranged parts of both lon- 
gitudinal and transverse correlation functions are very 
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This quantity is independent of the field, and in the 
NLsM approach is rather close to the constant valueiEl 
K = 1/2 given by bosonization for all ladder systems. 



VIII. COMPARISON WITH EXPERIMENT 

We turn to a comparison between the preceeding 
physical ideas and experiments performed on the two- 
chain, S — L/2 spin ladder material CuHpCl. The 
magnetizationtj illustrates the transition from the gapped 
regime of unbroken symmetry, and the nature of the 
broken-symmetry phase, while nuclear magnetic resp-i 
nance (NMR) measurements of spin relaxation ratesE2l 
reveal some characteristics of the spin correlation func- 
tions in both types of ground state. 

From Sec. VI, the magnetization takes the simple form 
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(8.1) 



where M (Eq. ( |6.lD ) becomes constant for the low tem- 
peratures of most interest. While the linear term is al- 
ways present, we have shown in Sees. IV and V that below 
a threshold field B* , where the system has a spin gap, it 
is canceled by the corresponding correlation-length term. 
As discussed in Sec. II, above an upper threshold (B c2 ) 
the magnetization will saturate at the value M s . Here 
we will for clarity treat saturation as a simple cutoff, 
rather than pursuing the more complex but fully system- 
atic method of Sec. II, which will be discussed elsewhere. 
Because the transition at B c2 is from an XY phase with 
strongly suppressed quantum fluctuations to a fully po- 
larized state with a spin- wave description, an exact treat- 
ment is of less interest in the current context. Finally, 



the 0(1/ N y S) term is included in Eq. (3T) as a reminder 
that large corrections are to be expected in any quanti- 
tative comparison of a "highly quantum" (small N y S) 
magnetic system with the semiclassical NLsM. 

Specializing to the parameters of CuHpCl, the ex- 
change constants deduced—from magnetization and sus- 
ceptibility measurementstJB are J'/ks — 13.2K and J/k B 
= 2.4K, whence J/[i B = 13.3T and J x /fi B = 3.6T. Tak- 
ing the simplest case of constant M , and the lower crit- 
ical field B c \ for onset where M{B C \) = 0, we obtain 
B c i = J jg[iB = 6.6T. The saturation field B c i is given 
from M(B c2 ) = M s as B c2 = (AS/g/j, B )J = 13.3T. These 
values agree well with a linear extrapolation of the mag- 
netization data at the lowest temperature in Ref. g, which 
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FIG. 6. Computed magnetization for the spin ladder sys- 
tem CuHpCl (see text). Circles are data from Ref. [5]. 



yields B cl = 6.8T and B c2 = 13. 7T. The computed mag- 
netization is shown in Fig. 6, where the dashed line in- 
dicates the validity limit (|6| > cir/a) of the calculation. 
The NLsM prediction of linear magnetization-aceears 
closer to the data than a repulsive boson modelEil — and 
other purely ID approaches, sueh as bosonizatior£3 and. 
simulations on a single laddcrE-tl (or a Haldane chain,E3 
which has the same universality class), all of which pre- 
dict souare-root cusps at B c \ and B c2 not evident in the 
data.H For a temperature T = 0.42K ~ 4^^ ^ s f° rm 
should be visible in the pure, Id systemEJ. This may 
be indicative of a weak interladder coupling, which also 
causes the real material to display 3D order at interme- 
diate B and low T. The gradient of the linear region can 
be written as 
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where Xaf = (ff/ i s) 2 /4>/a 2 is the Neel susceptibility per 
unit volume, the result for a classical antiferromagnet 
with effective exchange coupling J. The last term ex- 
presses the deviation from classical behavior as a func- 
tion of spin magnitude and the number of coupled chains. 
That the magnetization adopts such a quasi-classical 
form in the broken-symmetry regime, at least in systems 
with no topological term, may be ascribed in part to 
suppression of quantum fluctuations by the field (effec- 
tive beyond -B c i), and possibly also to the presence of 
higher-dimensional couplings. 

Further commentary is in order concerning the agree- 
ment between the NLsM results and experiment. The 
semiclassical model has recently been used with consid- 
erable success, also for the field response, oL&mall (spa- 
tially uniform) systems with rather small SS3 However, 
the apparent convergence of theory and experiment for 
the "maximally quantum" case of N y S = 1 in CuHpCl is 
in fact a somewhat fortuitous consequence of the param- 



eter ratio J/J±. The magnetization computed as above 
is valid for large B, and unlike the correlation length or 
gap is not sensitive to the coupling constants. Thus the 
model predictions for the isotropic ladder (J/ J' = 1), 
and for the Haldane chain (also N y S = 1), are little dif- 
ferent from the above, whereas the known results have 
much smaller gapped regions. Thus we stress again that 
the main value of the current approach lies-jji its quali- 
tative features. From this and other work,t3 it appears 
that N y S > 3 is sufficient for the NLsM to provide quan- 
titative accuracy, although the gap is then very small. 

The experimental comparison gives additional perspec- 
tive on the utility of the spin stiffness and 1/JV treat- 
ments (Sees. Ill and IV). We deduce from L* n that 
B* = 0.32J/g[iB — 4.2T, a value rather lower than B c \ 
above, and emphasize that in the regime between B* and 
B c \ neither treatment is reliable. Following the RG anal- 
ysis, the field scale 6* from the numerical solution gives 
B* = 1.6T, although prefactor corrections (above) can 
be expected. 

The other category of experiments performed on CuH- 
pCl is the measurement of NMR spin relaxation rates, 
which probe the spin-spin correlation functions discussed 
in Sec. VII. The experimental observable is the spin- 
lattice relaxation rate, given in d dimensions by 
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in which Fij are hyperfine coupling constants. This is 
obtained from the dynamic susceptibility 
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which in turn makes contact with the results of Sec. VII. 
By straightforward scaling, spin-spin correlation func- 
tions (T T Si(x, t)Sj(0, 0) ) decaying exponentially in 

1+ld space correspond to relaxation rates which are ex- 
ponentially activated in temperature, for T below the 
characteristic energy scale_A. This is observed in exper- 
iments for fields B < -B c i.lj Similarly, power-law spatial 
decays correspond to thermal power laws in relaxation 
rates: by performing separate Fourier transformations on 
x — ir and x + it and extracting the energy-dependence 
in the long- wavelength limit (k —> 0), if 
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whence from Eqs. (^2 7.3) the relaxation rate for the 
ladder system is 
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Ay and B± are constants related to the hyperfine cou- 
pling. An identical result was obtained in the Luttinger- 
Liquid formulation of Ref. p3. In addition to the acti- 
vated form in the gapped regimes at B < R-i and B > 
B C 2, the experimental measurements of 1/Tic3 show clear 
power-law divergence in the intermediate field regime, 
fully consistent with expectations (Sec. VII) based on a 
model of symmetry breaking by the effect of the field, 
on quantum fluctuations. However, because the dataE3 
cover at most one decade in temperature, while typically 
3-4 decades are required to identify the exponent within 
meaningful bounds, we may not attempt to make any 
quantitative deductions (Eq. (7.4)) concerning the cou- 
pling constant g. 



IX. SUMMARY 

In conclusion, the nonlinear a model description of 
quantum antiferromagnets in an external field provides 
a versatile framework for illustrating the breaking and 
restoration of symmetries. The case of 0(3) spins has a 
variety of physical realizations, which illustrate clearly 
the transition from the gapped, short-range ordered 
phase where full spin symmetry is maintained by quan- 
tum fluctuations, to a gapless, quasi-long-range ordered 
regime where the 0(3) symmetry is reduced to 0(2) by 
the field. 

The foregoing analysis is not restricted to CuHpCl, 
but applies to any system with a trivial topological term 
J2i-Pi (Eq- (2.2)). We comment here on its application 
to materials in two other classes. CFTDd is a planar 
S = | Heisenberg AF, of a type studied extensively in 
connection with the high-temperature superconductivity 
problem J3 and knowrfi to have AF order at T = 0. In an 
applied field one expects linear magnetization character- 
istics with the Neel susceptibility, beginning at B = 0, 
but not saturating because the superexchange interaction 
J = 9.4meV is beyond the reach of laboratory fields. 
Magnetization results on further recently-synthesized, 
planar Cu-eompounds with significantly smaller exchange 
constantsEH may show evidence of a logarithmic approach 
to saturation, as expected from spin-wave theory for a 
purely 2d system, and obtainable in the NLsM by inclu- 
sion of the total spin contraint as in Sec. II. The cases of 
most interest are those exhibiting the physics of the Hal- 
dane gap. The primary example is the S = 1 AF chain, 
of which NENP is considered a prototypical case (but for 
the complication of a large single-ion anisotropy). The 
present study yields the same qualitative features of a 
gapped regime with zero magnetization, followed by ap- 
proximately linear behavior towards saturation (which 
could not be achieved), as in experiment. a Quantitative 
accord with .experimental transition fields, and numerical 
simulations ,E3 does not match that of Sec. VIII for CuH- 
pCl. Finally, there has been considerable recent interest 
in the possibility of field-induced magnetization plateaus 



in certain systems: in a NLsM approach these may be 
expected, for example in S > 1 chains, when the field 
strength is such that the projected in-plane spin \Sn±\ is 
of integer amplitude, leading to a gapped phase. 
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APPENDIX A: RENORMALIZATION GROUP 
EQUATIONS 

The append ix presents the derivation of the integral 
forms (5.3,5.4) of the RG equations for the coupling con- 
stant and field. We begin with the Lagrangean for the 
NLsM in a magnetic field in the form of Eq. (5.1), 



Cl 



1 

2^ 



[(d^n) 2 -b 2 (l-n 2 z )+4ibn x n y 



(Al) 



where the parameters are defined in Sees. II and III, and 
the last term is obtained by partial integration of one 
component of the cross product. The staggered magneti- 
zation satisfies |n| 2 = 1 and (Sec. II) lies predominantly 
in the {x,y) plane, so that we may choose \n z \ -C 1 and 
specify the components of n by (n x = ^1 — n 2 z cos </>, 
n y = y^l — n 2 sm<j), n z ), with < cf> < 2tt. Using these 
variables, the cross product term gives n x h y — n y n x = 



(1 



and the gradient term becomes 
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(d,n) 2 = (d,n z ) 2 + (1 - nj)(d^) 2 + v ^ \' . (A2) 

i n z 

On making the substitution n z — y/g<7 z , where the cou- 
pling constant g is taken to be a small parameter (weak- 
coupling expansion), the Lagrangean ( |Al| ) is exactly rep- 
resented as 

C E = i(l - ga 2 z )[{d^) 2 + 2ib<j> ~ b 2 } + ±(d^) 2 
1 9 



2 ,(er z d M cr z ) 2 



2(l-.9^ 2 ) 

With the definition A' = (d^cj)) 2 
panded in the power series in g 



2il 



(A3) 
Ce may be ex- 



C E = &[A'- b 2 } + \{d^ z ) 2 - \(a z ) 2 \A' - b 2 } 



2.7 



-hiv^a,) 2 + y 2 a 2 z (a z d^a z ) 2 + 0(g 3 ). (A4) 



We proceed to one-loop order by neglecting those terms 
O(g) and higher, and by following the Wilson RG tech- 
nique of integrating out states of higher momenta around 
a finite upper cutoff.HLBlij In this process the renormal- 
izing effect of the fast Fourier modes of a z is taken into 



account, whereas the slow variable 4>(t, x) and its (small) 
derivative terms give no contribution at high momenta. 
We compute the quantity 



/= / Va z (k)e 



2 J d> 



2<l 



(A'-b 2 )+a z (-d 2 -A'+b 2 )a z +0(g) 



(A5) 



in the isotropic momentum shell 7A < |k| < A, where 
7 < 1 and A = n/a. Evaluation of the Gaussian integral 
yields 



I = 



n 



2~ 



k 2 - A' + b 2 



-±[A'-b 2 

e 2 9 



7A<|k|<A 



(A6) 



in which the small quantity A' permits the expansions 



Tr' \n{-dl -A'+ b 2 ) = Tr' ln(-^ + b 2 ) 

1 



and 



Tr' 



-d 2 u + b 2 ' 



:A' 



-Tr' 



Tr' 



-d 2 a + b 2 



.4' 



-d 2 + b 2 



A'. 



(A7) 
O(A') 2 

(A8) 



Sub stitu ting these results into the original Lagrangean 
Ce (A4) returns the new, effective Lagrangean on renor- 
malized lattice length scale a', 



£cff= H^ Tr '^W 



A' 



1 
2 
1 ( 1 



+ ^d,a z ) 2 --(a z ) 2 [A'-b 2 } 



2 Vffo 



b 2 -Ti'\n(-dl + b 2 ) 



(A9) 



Finally, by comparison of C c s with Ce , we may identify 
the renormalized coupling and magnetic field terms as 
expressed in Eqs. (5.3 pA). 



APPENDIX B: SPIN CORRELATION 
FUNCTIONS 

In this appendix we present the calculation of spin-spin 
correlation functions, in real space and inverse tempera- 
ture, for the gapless, high-field phase of the NLsM. The 
latter part of Sec. VII illustrates the transformation of 
these quantities for comparison with experiment. It is 
apparent from Eq. (fTTh that four such functions must 
be considered, corresponding to correlations of the stag- 
gered and uniform spin components in directions parallel 
(i,j = z) and perpendicular (i,j — ±) to the field. We 



reexpress the spin components n and 1 in terms of the 
angles ?9 and <f> (Sec. II), for which the propagators are 
known. 

The parallel response function for the staggered spin 
components is computed using 

(n,(x)n,(0)) = Jvnn z (x)n z (0)e- s W, (Bl) 

where n z = cos 9 = sin $ ~ $ for small angular deviations 
from the plane determined by the field. One expects the 
linear expansion to be valid for the lowest-order form 
of the correlation function, and may use it also in the 
approximation JI?(cos$) ~ / £*$■ Now 

Mx,tX(0,0)) = / 'l>0T?(x)0(O)exp [~±$ ■ G- 1 ^} 

= ] im o^4/ Wexp H"' G " 1 " +J '"}' (B2) 

in which ■& denotes a vector representing states at all 
x = (t,x), Gq 1 = —d 2 + b 2 is the matrix propagator 
for modes of 1?, and the current j(x) is taken to zero at 
the end of the calculation. Evaluation of the Gaussian 
integral yields 

A A 

lim Y~ Y~ ex P {l93- Goj} = 5G (x,0), (B3) 

J-+o 6j x djo 

from which the propagator component is given by reso- 
lution into plane wave states according to 



G (x,0) = (x| 



1 



10) 



1 



■£■ 



„ik-x 



-d 2 + b 2 ' ' LL T V k 2 + 6 2 



, (B4) 



where k = (w n , k). For the interesting limit of large dis- 
tance (by/x 2 + t 2 > 1), we conclude that the correlation 
function 



(n z (x,T)n z (0,0)) = 



A 



— b\/x 2j rT 



yby/x 1 + T 2 



(B5) 



where A is a constant, has the exponentially decaying 
form in space and inverse temperature expected from its 
origin in the gapped, out-of-plane excitations. 
i-The other component of the parallel response function 

i S y 

(*»(x)J.(0)) = (l?(*)I*(P)) 



AcLL T 



^e"' x (l-(n z (x)« z (0))). (B6) 



The latter term is a quantum correction arising from 
the constraint 1 • n = 0, which is unaffected by the 
field, and has the value (g/4c) (5 x ,o — (( n z(0)) 2 ))j «- e - 
is a local quantity which contributes only to a constant 
proportional to the square of the local magnetization. 
From Eq. ( |2.4| ), the expression for if may be simpli- 
fied by noting that the field term b = (0, 0, b z ) selects 



n z , which in turn introduces $, and thus only expo- 
nentially decaying contributions. The remaining term is 
(nAn) z = 0sin 2 i? ~ ^ only. With the exception of fur- 
ther constant terms ((n 2 (0))), which are neglected both 
here and in Eq. (7.2), we may restrict our considerations 
to 



(aS) 2 (l z (x)l z (0))=- 



1 



(4J) a 



{AJf 



(0(x)0(O)) 
^(0(x)0(O)>. (B7) 



We see now that we are dealing only with the gapless, in- 
plane modes, which can be expected to give a power-law 
correlation function at long distances, and the treatment 
is identical to the Luttingcr liquid. By the same tech- 
nique as above (Eqs. (B2 B^)), the correlation function 
for the azimuthal angles is given by 



(<Hx)0(O))=2. 9 BGo(x,O) 



(B8) 



with propagator Go = — d 2 for in-plane modes, and B 
a constant. Use of the result (expressed in rectangular 
coordinates with c retained explicitly) 



C ^-^ LU 2 


ikx+iujT ° ■*' 


T 2 


LL T ^ k 2 + (ujI 
c 3 


c) 2 " ~ 2tt [x 2 + c 2 
1 1 


T 2 } 2 


4-7T 


(x — icr) 2 (x + icr) 2 





(B9) 



leads directly to Eq. (|7.2| ). The prefactor is obtained 
from the identity c/2J = a/ g, while further identifica- 
tion of 7r/2g with K ( [7.4| ) yields the result of Ref. p^ . 
Because we are not interested in finite-size properties for 
discussing the asymptotic behavior of c orre lation func- 
tions, we ignore the zero modes (Eq. ( |B9| )) occurring 
in the Luttinger liquid. Their presence is connected with 
decompactification of <f> and the associated winding num- 
ber (Sec. VI), and their treatment discussed in Ref. EJ3. 
Moving to the correlation functions normal to the ap- 
plied field, that for the staggered spin components has 
the simple form 

S = (n+(x)n_(0)) = (exp(z0(x) - i<p(0))). (BIO) 

Defining the function j(x') = <5(x — x') — #(x'), and us- 
ing the equalities Go(x, x) = Go(0, 0) and Go(x, 0) = 
Go(0,x), permits one to write 

3 = J P0ex P {-^0 • Gq V + ij ■ 4>) 
= (DctGo 1 )- 1 / 2 exp{-igj.G j} 
= G'exp{- g [G (0,0)-G (x,0)]}, (Bll) 

where C is a constant. The propagators are evaluated 
as above, to give 



6 °(0,0) = ^£p = ^(£] 1B12, 



and 



G (x,0) 



I ^ e jkx 1 , {(2it) 2 (x 2 + t 2 ) 



E 



LL T ^ k 2 4tt 

k#0 



In 



C 2 



(B13) 



both of which are cut off by the system length scale 
(Sec. V) C = min[L,L T ]. Finally, 

_-2_ _A. 

~ / C\ 2tt / t 2 + T 2 \ i-JT 

(n+(x)n_(0)> = C^J \~^-) > ( B14 ) 

for altered constant G, returns the result in Eq. (|7.3|). 
That this is the only power-law contribution to 
(S+(x)S-(0)) is easily seen from 

(i£(x)£(0)) = (e** (x) -^ (0) 1?(x)i?(0)} 

= -(e^ (x) - i< * (0) )[^( t ?(x)^(0))], (B15) 

whose dependence on •d guarantees an exponential decay 
of correlations. 
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